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We study the effects of time-independent nonequilibrium drive on an open 2D electron gas sys- 
tem coupled to 2D longitudinal acoustic phonons using the Keldysh path integral method. The 
layer electron-phonon system is defined at the two-dimensional interface between a pair of three- 
dimensional Fermi liquid leads, which act both as a particle pump and an infinite bath. The 
nonequilibrium steady state is achieved in the layer by assuming the leads to be thermally equili- 
brated at two different chemical potentials. This subjects the layer to an out-of-plane voltage V 
and drives a steady-state charge current perpendicular to the system. We compute the effects of 
small voltages (V <C lod) on the in-plane electron-phonon scattering rate and the electron effective 
mass at zero temperature. We also find that the obtained nonequilibrium modification to the acous- 
tic phonon velocity and the Thomas-Fermi screening length reveal the possibility of tuning these 
quantities with the external voltage. 
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I. INTRODUCTION 

When considering steady states in a non-adiabatic 
closed driven system it is necessary for one to specify 
dissipation effects. Intuitively, bulk heating by an ex- 
ternal drive must be compensated by bulk dissipation or 
some form of thermal contact with an infinite bath so as 
to prevent gradual heating. In the absence of a clear dis- 
sipative mechanism within the system one must consider 
some form of coupling between the system and an infinite 
reservoir through which heat can escape. However, it is 
difficult to model an infinite bath that would dissipate 
heat at every location in a 3D bulk system. Recent the- 
oretical worksi 2 -^^ in electrically driven steady-state 
systems investigated universal scaling behaviour in trans- 
port quantities near various quantum critical points. In 
the work by Dalidovich et a&, a heat sink was consid- 
ered by adding a phenomenological dissipation term in 
the Lagrangian as suggested by Caldeira and Leggett- 
In the work by Green et a&, a heat sink for the itinerant 
electrons was assumed to be provided by the underlying 
lattice. 

In the treatment of these closed nonequilibrium sys- 
tems, it is often the case that precise details of the cou- 
pling between such systems and their environment are 
not known, and one is often reduced to describing these 
effects phenomenologically in an ad hoc fashion. An al- 
ternative method of treating the heating problem in a 
system is to begin with a theory which explicitly includes 
couplings between the system and its external reservoirs. 
In an open driven system these reservoirs naturally act 
both as a source of nonequilibrium drive (particle pump) 
and a heat sink. A vast number of theoretical works 
on open driven systems have been conducted in meso- 
scopic physics. Examples can be found in quantum dot 
system o 10 : 11 ' 12 ! 13 where leads that couple to the dot can 
be envisaged as the reservoirs. 

Though the role of these reservoirs is vital, it is de- 
sirable to have a final effective theory formulated by de- 



grees of freedom associated with the central active sys- 
tem. From a mathematical viewpoint, this requires one 
to trace out reservoir degrees of freedom from the start- 
ing theory. This process of tracing out bath degrees of 
freedom is most straightforwardly done in the language 
of functional integrals. In the context of nonequilib- 
rium systems the Keldysh path integral method is a suit- 
able technique-'-' 9 . The formalism has been used pre- 
viously in characterizing zero-dimensional open systems, 
like quantum dots, subject to charge current o 11 ' 13 ' 14 . In 
the context of extended open systems the formalism has 
recently been applied to a two-dimensional itinerant elec- 
tron system^ and a microcavity polariton syste m 16 ' 17 . 

In this work we consider a steady-state open 2D metal- 
lic system which is electrically driven. We model the cen- 
tral metal layer in a jcllium model in which conduction 
electrons couple to 2D longitudinal acoustic phonons. 
Two metallic leads, which sandwich the metal layer 
(FigfTJ) , are in thermal equilibrium at two different chem- 
ical potentials, and establish an out-of-plane charge cur- 
rent through the layer. While the electrons are driven out 
of equilibrium by its direct coupling to the charge current 
the phonons are out of equilibrium due to their coupling 
to the electrons. Using the Keldysh formalism we inves- 
tigate the effects of nonequilibrium perturbation on the 
properties of electrons, phonons, and their interactions. 

We now summarize our main findings. First, we find 
that the zero-temperature phonon velocity is modified in 
the presence of a nonequilibrium voltage, 

Sc(V) ~ (AT)V. (1) 

V = (1l — is the nonequilibrium voltage defined as 
the difference between the two chemical potentials of the 
leads, and Ar = Tl — describes the asymmetry in 
the coupling strengths of the central metal layer to the 
two leads. We note that the sound velocity can be tuned 
by varying voltage V, and may be increased or decreased 
according to the polarity of the voltage. Screening prop- 
erties of electrons are also modified by the voltage. In 
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particular, we note the nonequilibrium correction to the 
Thomas-Fermi wavevector, 

Sk s ~ (AT)V, (2) 

which shows how the voltage can modify the distance 
scale beyond which the Coulombic disturbance of the ions 
is effectively screened by the conduction electrons. We 
also consider modifications to in-plane transport in the 
presence of out-of-plane current. In particular, the out- 
of-equilibrium electron-phonon scattering rate at zero 
temperature, for voltages much smaller than the Debye 
frequency (V •< wp), scales with voltage as 

— V 3 . (3) 

T~el— ph 

We note that in thermal equilibrium, the scattering rate 
scales with respect to temperature with the same power 
law. The reason for this, as discussed in the main parts 
of this paper, can be found in the similar responses of the 
in- and out-of-equilibrium electron distribution functions 
to temperature and voltage, respectively. Corrections to 
the electron mass enhancement factor, A = with 
respect to voltage at zero temperature is 

SX(V) ~ V\ (4) 

with the finite-temperature correction scaling with the 
same exponent. Eqf4] implies that voltage can be used to 
tune the effective mass of the conduction electrons. 

The paper is organized as follows. We present and 
describe our model in Sec lII Al and derive an effective 
phonon action using the Keldysh formalism in Sec lII Bl 
In Sec lIIII we compute effects of nonequilibrium pertur- 
bation on the phonon velocity, electron-phonon scatter- 
ing rate, and electron mass enhancement. We conclude 
in Sec lIVI Details of the calculations in Sec lII Bl are pro- 
vided in Appendix O 



II. THEORY 
A. Model 

The geometry of our system consists of two 3D elec- 
tron reservoirs, or leads, which are in contact at a two- 
dimensional interface. We define this two-dimensional 
interface as our central metal layer of interest (FigfTJ. 
We split the Hamiltonian into three pieces: H = Hr, + 
Ht + Hi ayer , where Hi and Hi ayer describe the leads 
and the metal layer respectively, and H t models the tun- 
neling between the leads and the layer. We view the 
leads as noninteracting 3D electron gases. They are both 
in thermal equilibrium at temperature T but may pos- 
sess different chemical potentials fii and [1r, thus driving 
a steady-state current in the direction perpendicular to 
the metal layer. A continuum of states is assumed in 
the leads, occupied according to the Fermi distribution 



left lead /i L 2D metal 




FIG. 1: A schematic representation of the system. A 2D 
metal layer is sandwiched between two 3D leads. The leads 
are assumed to be noninteracting electron gases in thermal 
equilibrium at temperature T but may be held at two different 
chemical potentials. 



f°(x) = [l + e"' 1 -"-'] \ where a G (L,R) labels the 
leads. Hl is then given by 

Hl= £k ll> fc ^ C k||,fc J! , ( T,Q C k||,fc^o',Q- (5) 

cr,a = L,R 

The z-axis of our coordinate system coincides with the 
normal of the metal layer, which is defined on the z = 
plane. k|| denotes the component of a momentum vector 
parallel to the z = plane and k z is its component normal 
to the plane. Hereafter, we drop the "||" label on all 
momenta for brevity, and we take Ti = 1. 

The form of our tunneling Hamiltonian H t assumes 
conservation of spin a and parallel momentum vector 
k of the tunneling electrons, thus making a and k 
good quantum numbers throughout the entire lead-layer- 
lead structure. For computational simplicity, we assume 
energy-independent tunneling matrix elements but main- 
tain their lead-dependences in order to describe possible 
asymmetries in the lead-layer couplings. These assump- 
tions can be summarized in the following form for the 
tunneling Hamiltonian, 

<r,a = L 1 R 

Finally, the central layer is a 2D metal whose Hamil- 
tonian is, 

Hiayer = ^ e ^A.,a d K" + f^k + \ j 

k,cr k ^ ' 

k,q, a 
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V£ e = 2Tre 2 /q is the 2D Coulomb potential. Our theory 
assumes a jellium model for the electron-ion interaction 
and consequently only describes a coupling between lon- 
gitudinal acoustic phonons and electrons. The effects 
of transverse acoustic phonons are not included in this 
work because most important properties of a metal can 
be understood by observing the coupling of electrons to 
changes in the background ionic charge density which are 
described by acoustic phonons. The discussion of optical 
phonons may be neglected given that our work only con- 
siders monatomic lattices. In the 2D jellium model of the 
electron-phonon system, the unscreened electron-phonon 
coupling can be shown to be 



5q 



2Mw. 



(8) 



where e q is the polarization vector of the longitudi- 



nal phonons for each q and V% 



ZV* 



u) p (q) 



yj2imiZ 2 e 2 q/M is the 2D ionic plasma frequency and 
u q = 6 q + 6^ q is the dimensionless phonon displacement 
operator. M, Z, and rij denote the mass, valence number, 
and number density per unit area of the ions respectively. 



B. Effective Nonequilibrium Phonon Action 

In this section, we obtain a description of the lattice 
vibrations in the central metal layer. The metal layer 
consists of two main components: interacting electrons 
and phonons. When the layer is coupled to two elec- 
tron reservoirs with different chemical potentials, one can 
think of the bath with the higher (lower) chemical poten- 
tial as a particle pump (sink) to the central active layer. 
Clearly, if the rate of thermalization within the layer is 
much smaller than the rate of tunneling between the layer 
and the leads, the electrons in the layer ceases to follow 
the Fermi-Dirac distribution. This nonequilibrium na- 
ture of layer electrons also influences the properties of the 
phonons due to electron-ion coupling. Therefore, the sys- 
tem must be described using a nonequilibrium treatment 
of the layer electrons, phonons and their interactions. We 
use the Keldysh formalism to provide a unified descrip- 
tion of the two-dimensional electron-phonon system both 
in and out of equilibrium. 

The Keldysh formalism requires each field to possess 
two values at each point in time, one on the forward 
branch and another on the backward branch. As a con- 
sequence, electron (G) and phonon (D) Green functions, 
along with their respective self-energies (E,H), must be 
represented by matrices 



G 
t 



G R G K 
G A 

E A 



,D = 



,n 



D A 

o n A 



(9) 



For the many-particle Hamiltonian introduced in 
Sec lII A[ the Keldysh generating functional, Zk, can 



be obtained by a familiar real-time coherent state func- 
tional integral technique with the only modification aris- 
ing from performing the real-time integral over a time- 
loop contour^. After integrating out the lead degrees of 
freedom and introducing a Hubbard-Stratonovic field 4>, 
which decouples the two-body electron interaction term, 
the effective Keldysh generating functional yields 



T>[d, d, u, u*, <j>, ( / ) *] e istf f Mu, v ;w] 
V[d,d,u,u*, <M*] ( 10 ) 



x e 



i^S'^ld.Au.u* ,4>,<p*] + S'j£[<p,<p*] + S p K h [«,«*]) 



where 

iS e ^[d,d,u,u*,(f>, (f>*} 



2i(dtY J (4(k',i) #(k',i)) 

J k,k'ff 



G- 1 (k',k,i) + <l>(k'-k,t)| ( "f 



4(M) 

d»(M);> 



2i J dtJ2(r l *(Kt) <F*(M) ) 



o £ 
& 



d (M) 
^(M) 



(ii) 



(12) 



2i I dt^2(u cl *(k,t) u'*(k,i) 

« 1- 



(13) 



In the absence of leads, hence for a closed system in equi- 
librium, the inverse Green function matrix, 



G^(M)<5 k , k [G- 1 (k,t)] K 5 k , lt . 
G^(k,t)5 k > k 



(14) 

describing the propagation of electrons in the central 
layer has the usual free-fermion form. However, these 
Green functions are modified due to tunneling between 
the layer and the leads. For energy-independent tunnel- 
ing amplitudes tL,R, the dressed Fermion Green functions 
were calculated elsewher o 13 ' 14 , 



G R (K) 
G K (K) 



1 



[G A (K)]* 



UJ — £ k + iT 

-2 l E Q r Q tanh(^ 



(15) 



2T 



(^-e k ) 2 + r 2 



Here, r = Tl + Fr, and r Q = nvt^ is the rate at 
which electrons escape from the active layer into lead a. 
v = m/ir is the 2D density of states (including both spin 
projections) for electrons at Fermi energy. We have also 
introduced an abbreviation, K = (k, to) . The $ matrix 



in EqQj] encodes the coupling of layer electrons to both 
the Hubbard-Stratonovic field and phonons. It is given 
by 



*(M) 



$ c/ (k,i) $ 9 (k,t) 
$«(M) $ ci (k,i) 



where 



$ cZ (k,i) =e^(k,t)- 3kW d (k,t), 
$«(k,t)=e^(k,t)- 9k u«(k,t). 



(16) 



(17) 



The inverse phonon Green function matrix, in the ab- 
sence of electron-ion coupling, is given by 







[fl» B (M) 



[r»^(k,t)]-i_ 

^(M) 



(18) 



Z?o (k, t) describes the collective modes of the lattice ions 
as they oscillate in a negatively charged static back- 
ground of electrons at the 2D ionic plasma frequency, 
u) p (q). The bare retarded and Keldysh propagators have 
the usual form, 



D o(Q) = 



2uj p (q) 



D o(Q) = -27ricoth 



Q — u)p{q) + isgn(fl)5 



2T 



(19) 



X [8(il + u p (q)) + S(il - u p {q))} , 



where Q — (q, f2) is a phonon momentum-energy 3- 
vector. 

The next step is to integrate out the layer electron 
fields. We use the standard RPA treatment for both 
Coulomb and electron-phonon interactions. A schematic 
representation of the resulting action is 

iSx f [uci, u q , (f)cl,(/>q\ 

= 2iJ-^J^ T *(Q)D* 1 (QWQ), (20) 

where 



$ T (Q) = ( <MQ) *,(Q)) 



(21) 



The detailed form of this action can be found in the Ap- 
pendix [AJ The propagator matrix D^(Q) represents the 
effective electron-electron interaction. Two electrons can 
be scattered by the unscreened electron-electron interac- 
tion (wavy line) or by sending a phonon from one to the 
other (dashed line) (see FigEJ) . In the RPA approxima- 
tion, we consider diagrams where one or more polariza- 
tion bubbles are inserted in a chain. Both the interaction 
to the right and left of each of the bubbles may either be a 
Coulomb line or a phonon line. All of these diagrams are 
summed exactly in the path integral formalism by solv- 
ing for the inverse of the propagator. Within the RPA, 




FIG. 2: A diagrammatic representation of the propagator 
The propagator, which corresponds to the curly line in the fig- 
ure, describes the effective interaction between two electrons. 



the effective phonon action can be obtained exactly by 
integrating out the Hubbard-Stratonovic field <f>: 



iS e ^[u c i,u q ] 



2i 



d 2 q f <Kl _, T 



(2nf 



2vr 



u T *(Q)D- 1 (Q)u(Q), (22) 



where 



u T (Q) = (u C i(Q) u q {Q) 



(23) 



As shown in the Appendix [Al the inverse phonon propa- 
gator matrix Z) _1 possesses the bosonic causality struc- 
ture in Keldysh space. Its components are given by 



D- R 1 (Q,V) = 2 



[b-\Q,v)] K 



= 2 



n ° R[Q) e(Q,V) 



-u^k \ 9ci \ 2 Ti K (Q,V) 



(24) 



i D o (Q)r - 



e(Q,V) 



(25) 



where the nonequilibrium polarization functions are 
given by 



n R (K, V) 



[G K (K + Q)G A {Q) 



+G R (K + Q)G K {Q)} = Ii A *{K, V) (26) 



U K (K,V) = ^f [G K (K + Q)G K (Q) 

+G R (K + Q)G A {Q) + G A (K + Q)G R (Q)] . (27) 

Here, f Q =J d 2 qdQ./ (2ir) 3 . Voltage V is defined by V = 

Mi — l l R- e (QiV) denotes the nonequilibrium dielectric 
function which is defined by 



e(Q,V) = l-V^ e U R (Q,V). 



(28) 



We see in Eqs l24l25l that the voltage dependence of 
nonequilibrium phonon propagators is encoded in the 
self-energy part; this is consistent with the fact that 
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phonons are driven out of equilibrium only through its 
coupling to the electrons. The Keldysh phonon propa- 
gators (Eqs l24l25|) and the nonequilibrium polarization 
bubbles (Eqs l26l2"Tj) will now be used to calculate var- 
ious properties of our nonequilibrium electron-phonon 
system. In particular, we compute the effects of voltage 
on the phonon velocity, electron-phonon scattering rate 
and electron mass enhancement, and compare the results 
with the corresponding results in thermal equilibrium. 



III. PROPERTIES OF THE NONEQUILIBRIUM 
ELECTRON-PHONON SYSTEM 



A. Preliminary Considerations 

There are four energy scales in our system: temper- 
ature/voltage T/V, the Debye frequency ujd, the layer- 
lead scattering rate T and the average Fermi energy of 
the two leads \i. Clearly, the mass difference between ions 
and electrons lead to a large discrepancy between the dy- 
namic time scales of these species, i.e. un -C u. Due to 
this difference, we are justified to employ the adiabatic 
approximation for the electrons in which we assume that 
electrons take on the configurations they would have if 
the ions were frozen into their instantaneous positions. 
In this approximation, a frequency-independent (static) 
dielectric function can be used. 

We assume that the layer-lead scattering rate satisfies 
wd C T C /i. Working in this limit validates an ex- 
pansion with respect to parameter, Q/T <C 1, where fi is 
the typical energy scale of the excitations. The limit is 
consistent with the necessity of a strong layer-lead cou- 
pling in order for electrons to remain out of equilibrium 
midst other thermalizing effects within the layer. In all 
of our calculations, we assume max{T, V} to be small 
compared to all other energy scales, or more precisely, 
max{T, V} <C Wd . Another observation one can make 
is our use of a jellium model for the ions. This sets a 
natural cutoff in momentum at 1/a ~ qf, where a is 
the lattice constant, because any details corresponding 
to length scale a or less is washed away by this contin- 
uum treatment of the lattice. 

In the long-wavelength (g <C kp), small-frequency 
(f2 <C T) limit, and for small voltages (V <C T), the 
polarization functions Eqs l26l2"71 are given by 



n R (q,n) « - Cl (v) 



■ [£) -fcfcWp (29) 



r a r/3 



r 2 



x>-,,ii i r n + ^"^) n + ^" M/? . 00) 

Here, /i = (/jl + /Ur) /2 is the average chemical potential 
of the two leads. All coefficients (ci(V), C2{V), cs(V), 



and C4) are positive, real, and satisfies c±(V) ~ Cs(V) ~ 
C4 ~ v 3> C2(V). D(fj) is the density of states at /x in the 
presence of leads, and is given by 



D{ji) = - 

7T 



7r , fa 

2 +tan (r 



(31) 



B. Nonequilibrium Bohm-Staver Relation 

An explicit expression for the dressed retarded phonon 
propagator calculated in Eqf5U is 



D R (Q) = 



Q 2 



2cj p \g^\ 2 n R (Q,V) ' 



(32) 



which in the long-wavelength, static limit is reduced to 



D R {Q) 



Q-2 



2ire'm R (Q=0,V) 



(33) 



The pole of this propagator gives the dispersion of the 
non-attenuating longitudinal acoustic phonon mode, 

fi* = c 2 (V)q 2 , (34) 

where the voltage-dependent phonon velocity is, 

Z 2 n,. 



c\V) 



MILr(0,0,VY 



(35) 



At zero temperature, the static, long-wavelength polar- 
ization function as a function of voltage is 



a 

v 1 Ar ^ 



-DO*) 



^D(fi) 1 + ^ 



r r 



o 



V 2 



(36) 



where Ar = Tl — Tr. From Eqs l35|36"l the nonequilib- 
rium phonon velocity is given by 



c 2 (V) 



1 Ar v 



2 7 rZ)( At )i + ^ r r 



where 



4 



Z 2 n t 
AID {n) 



(37) 



(38) 



is the equilibrium phonon velocity as obtained by Bohm 
and Staver— . We notice in Eq[37] that the acoustic 
phonon velocity can be tuned by voltage. In particu- 
lar, depending on the tunneling asymmetry and polarity, 
the change in the velocity with respect to its equilibrium 
value may be positive or negative. 

The equilibrium phonon velocity Eq |38l may not be 
in the most familiar form. A more common form of 
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the equilibrium phonon velocity is derived from requir- 
ing charge neutrality n e — Zrii and the relationship of 
ratio n e /D(fi) to the Fermi energy. In the considered 
nonequilibrium system such relations do not hold due 
to charge accumulation and an absence of a well-defined 
Fermi energy in the metal layer. However, Eqs l35l38l sug- 
gest that both equilibrium and nonequilibrium phonon 
velocities have the form of Eq l38l where nonequilibrium 
effects are encoded in the generalized density of states, 
—Hji(O,0,V). In equilibrium, the density of states is 
evaluated at one single chemical potential fi which is lo- 
cated at the Fermi distribution edge. At finite voltage, 
the electron distribution function can be shown to have 
a split-step shape 

T e V) = z/V - Ml) + (1 - x)f(u - fjL R ), (39) 
where f° denotes the Dirac-Fermi distribution and x = 



FIG. 3: Nonequilibrium steady-state electron distribution 
function for our system with fiL — = V > and x — Yl/T 
is shown above in bold (in the absence of electron-phonon 
interaction effects). The Dirac-Fermi distribution function is 
represented with a thin line. 



/r. The nonequilibrium analog of the density of states 
at Fermi energy is the sum of density of states at [II 
and hr weighted by the respective lead-layer coupling 
strengths, and is given by: 



n fl (0, 0, V) = xD{ii L ) + (1 - x)D{jjl r ) 



(40) 



This voltage dependence in — IIr(0, 0,V) is the origin 
of the voltage dependence seen in the nonequilibrium 
phonon velocity. 



l/r e i- P h ~ T 2 . This is what is expected in a closed sys- 
tem. In the considered open system, electrons can scatter 
into one of the leads frequently so that the chance of an 
electron encountering a phonon and scattering with it 
is decreased. Therefore, we expect a further decline in 
the electron-phonon scattering rate. We then extend our 
calculations to the zero-temperature finite voltage case 
where the scaling behaviour of the scattering rate with 
respect to voltage is determined. It is crucial to note that 
the behaviour of phonons depends most heavily on the 
nonequilibrium distribution function of the electrons. In 
this respect, the T — 0, V ^ situation is similar to 
the T 7^ 0, V = situation in that the latter involves a 
smearing of the step distribution while the former results 
in a split-step distribution as shown in Figj3] There- 
fore, phonons are subject to electrons whose distribution 
behaves similarly in both situations, and one may in- 
tuitively predict the scaling behaviours of the electron- 
phonon scattering rate with respect to T and V to be 
similar. We now quantify this possible similarity. 

The starting point of this calculation is the one-loop 
self energy of the electron (FigHJ). We may invoke the 




k, u k + q, uj + f2 
FIG. 4: One-phonon self energy of the electron. 



Migdal theorem and ignore vertex corrections. Formally, 
the self-energy in Fig 0] is given by 



Z R (K) = 2i 



lff(q) 



HQ>V)\ 



2 (G (K + Q)D (Q) 
+G K (K + Q)D R (Q)) . (41) 



In the static, long- wavelength (q <C fop) limit the dielec- 
tric function is purely real and is given by 



2TTe 2 Cl {V) _ k s (V) 



(42) 



C. DC Electron-Phonon Scattering Rate 

Here, we discuss the electron scattering rate due to 
longitudinal acoustic phonons. In equilibrium, at tem- 
peratures well below the Debye temperature, phonon 
wavevectors are of order T/cq or less. Then, within the 
wavevector space of phonons that are permitted to be ab- 
sorbed or emitted by conservation laws, only a subspace 
of linear dimensions proportional to T can actually par- 
ticipate. In two spatial dimensions, this subspace is one 
dimensional so the phase space of allowed wavevectors is 
proportional to T. Combining this with the fact that the 
screened electron-phonon coupling |g sc | 2 ~ q ~ T, we 
expect the electron-phonon scattering rate to decline as 



where ks(V) = 27re 2 Ci(V) is the 2D nonequilibrium ana- 
log of the Thomas-Fermi wavevector. At V — 0, c\ re- 
duces to D(/j,) and the 2D equilibrium expression for the 
Thomas- Fermi wavevector, kg = 27re 2 D(/i), is recovered. 
We notice that the voltage dependence observed in the 
Thomas- Fermi wavevector indicates that voltage can con- 
trol the effectiveness of conduction electrons in screening 
Coulomb interactions between ions. Using H K in the 
static approximation, 



n 



~ ic ^ coth i 

a, /3 



2T 



r 



(43) 
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the dressed phonon Green functions in the static limit and 
can be computed as follows: 

D a (Q) = ^ (44) 

I 



D K (Q) = 



\D-\Q)] 



K 



{D*iQ)}-i[DA(Q)Y 



n 



2*coth( ^)l m [D^(Q)}- 1 



9ci \ 2 Tl K (Q,V) ImD R (Q) 
\e(Q,V)\ 2 Jm[D*(Q)]-i 



, 

-TTl COth — - 



2TJ Lo q {V) 



(6(n-u g (V))-6(n + u g (V))) 



+ 



-2ic 4 u p uj q (V) |.g q | Y. a . 



8 r 2 



£ coth r \ 



2T 



}(V) \Vl 2 - w 9 (y) 2 + «s 5 n(ft)<5| 



(45) 



In equilibrium, II = 0. In this case, by comparing the 
first term of Eq05] with EqfJ31 we can easily verify that 
the fluctuation dissipation theorem is satisfied. 

Now one can use Eq[4T] and the relevant Green func- 
tions to calculate the general DC scattering rate valid 
both in and out of equilibrium. Defining the static 
screened electron-phonon coupling constant as 

\g sc HV)\ 2 ^ 2 ^ q f7 T/V (46) 
e 2 q {V)uj q (y) 

I 



we obtain 



£ fl (k = k F , oj = n, T, V) = Z R (T, V) 

= Z R (T,V)+H R (T,V), (47) 

where 



1 



fl — vfQ cos 9 + iT 



-iri coth 



2T 



(5(n-Lj q (v))-s(n + u q (v))) 



+ 



-2ic^ v Lo q (y) | ffq f £ Q>/3 %^ coth 



^ i1 + tJ, a -flf3 \ 



e 2 (V) \fl 2 - uj q {V) 2 + isgn(fl)5\' 



(48) 



and 



E?(T, V) = 2i 



d 2 q dfl 

(2tt) 2 2^ 



5sc(q, V)r 



-2i£ Q r a tanh(^^) 



(Q - v F qcosdf + T 2 ^ 2 - w g (F) 2 + isgn{Q)S 



(49) 



Here, we have linearized the Fermion dispersion and Vl.r = ±V where V = [II — fJ-R > is assumed. Eqs l48l49l are 
the central results in the remaining sections of our work. In equilibrium (i.e. (J.l = (J-r), Eqs l48l4"9l reduce to 



Sf(T) 



qdqdB ,^,2 c «th (^) [S(fl - u q ) - 5(0, + w q ) 



(2tt) : 



-dQ \g sc (q)\ 



(fl — VFqcos9 + iT) 



(50) 



and 



Sf (T) = 4r 



qdqdO dfl . l2 
72^2^ c(q)l 



tanh (^) 



{fl-v F q cos 9) 2 + T 2 ^ 2 - to 2 + isgn{fl)5' 



(51) 
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In the presence of voltage (at zero temperature) , 

(27T) 



1 



£1 — vpq cos + iT 



- [ir S gn{n) {S(Q - u q (V)) - 8 (to + u q {V))) 



2c 4 w p uj q (V) |g q | 2 (sgn(Q + VQ - sgnjfl - V)) 

&(y) \n 2 - iv q (V) 2 + isgn(n)5\ 2 



and 



Xf(V) = 4 



d 2 q <m 
(2n) 2 2^ 



ff sc (q, V)| 



£ a r aSff n (n - v a /2) 



(Q - v F qcos6) 2 + T 2 ^ 2 - u q {V) 2 + isgn(Q)5' 



(52) 



(53) 



1. Thermal Equilibrium: T 0, V = 

In equilibrium scattering rate calculations, tempera- 
ture T defines a natural cutoff for the momentum inte- 
grals. This is because only phonons of frequencies com- 
parable or less than T can be absorbed or emitted. For 
T <C u>Di the phonon dispersion obeys tu q = c^q which 
implies that phonon wavevectors are of order T/cq or 
less. We now proceed with the integrals in Eqs l5QI5H 
The angular integration can be done straightforwardly. 
For atj£ < 1, 



(16 



1 



o 2tt lu + fi — v F q cos 6 + iT 



-i/o(«) + /i(«)- 



n 



where 



fob) 



1/2 ' 



(54) 



(55) 



conducted, the dielectric function and fo can be approx- 
imated by 



1.0 1 

<9) « ~f Ms) « ~- 



(59) 



Similarly, we get 



^/^tanhf^ 
Teg 7 2tt 1 2 



7 2 



(60) 



where 7 2 ~ O(l). Therefore, in thermal equilibrium, 



1 _ G 7 e « 

T el-ph 1 c 



T 1 



(61) 



where j eq ~ 0(1) 



1 



v r 



Using this result one obtains, 

rT/c 



3/2 ' 



(56) 



(I) 



GT 3 n f 1 y 2 dy 

■2 / ^— ^coth 



Teg "7o 2tt 



7i 



Tel 



T , 



where ~ 0(1) and 



G(V) 



\9sc(q,V)\ 



(57) 



(58) 



is a real quantity independent of Q. We have also used 
the fact that in the region of q over which the integral is 



2. Out of Equilibrium: V / 0, T = 

In the noncquilibrium case, the momentum cutoff is 
set by the voltage, i.e. q c ~ V/c(V). Therefore, the DC 
scattering rate is at least of order V, and thus, we may 
set V = for the quantities in the integrand such as the 
phonon velocity and dielectric function. This should al- 
low us to obtain the scattering rate to leading order in 
V. When the first integral Eqj52]is naively evaluated we 
find that the second term, which becomes finite out of 
equilibrium, gives a divergent contribution to the scat- 
tering rate. This divergence can be cured by dropping 
the static limit assumption and expanding the dielectric 
function to lowest order in frequency, 



e(Q) wl + d 
Recall that for \i ^> T, 



ic 3 - 



2ire 2 ft 



Cl 



c 4 ps 2c 3 



(62) 



(63) 
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Therefore, we get gJf _ 2 |g q | 2 H K (Q, V) tj 

(64) 1 ^ j 

where fcjy = 27re 2 ^. This new expression for the dielectric 
function assumes that typical excitation energy scale is 
smaller than T. The phonon dispersion then becomes 
dissipative: 

LQ 2 (q) f n 

g e(Q, V) q \ 2T " 



where oj q is phonon dispersion in the static limit. The 
modified phonon Green functions are, 

D R (Q) = - (66) Now we are set to compute the first self-energy term, 

ft 2 - uj 2 {V) +iuj 2 q {V) w 

I 



JmEf (V) = 4 



qdq dfl 
~2n 2tt 



ff S c(q)| 4 /o(<?)- 



-ImIL K (Q) 



(1 



n 2 

4T 5 



1 - 



n 2 



+ 



n 2 
W 1 



(68) 



We see that majority of the contribution in the integral 
comes from fl ~ u) q . So one may cutoff the Sl-integral by 
V < T. If we cutoff q at q c ~ f , f2 c - V < I\ The 
static expression for the Keldysh polarization function 
must be made consistent with the approximate dynamic 
dielectric function. For < V, 



-Imn K (Q)^D 



n 



r 3 



r 2 



2T L T R 

r 3 



V 



After a series of integrals, we get 

c 2 D 01 r q 



- ImSf(V) = 



V 2 r A 
1 c 



1 



2r L rj 
r 2 



V' 4 



(69) 



(70) 



where 7™ e<? ~ 0(1). This gives us a fourth order contri- 
bution and will become a subdominant correction to the 
second integral as we will now show. The second integral 
from Eg 1531 can be done straightforwardly: 

1 r V/co qdq 

{sgn(u q - V/2) 
GV 3 f 1 y 2 dy 

'1/2 



ImZ R (V) 



\9sM)\ z T 
sgn{oj q 



V/2)) 



1 Cq 



2tt 



7 2 



G 7 2 



neq 



Tel 



-V A 



(71) 



where eq 
we obtain, 



O(l) 



Therefore, to lowest order in voltage 
G 7 " e ' - 



1 



neq 
el—ph 



rc 3 



(72) 



where y neq ~ 0(1). 

We can summarize the low-T/low-V electron-phonon 
scattering rate both in and out of equilibrium: 



1 



7~el—ph 



1 Cq 



Where \ X=\V\ : T = 0,y^0 ' (73) 
where ■yx is some real constant of O(l). 



D. Electron Mass Enhancement 

We now compute the effective on-shell mass correction 
and investigate the effects of temperature and voltage to 
lowest order. 



1. Thermal Equilibrium T / 0, V — 

From Eqs l5QI5Tl we obtain, 
i?e£f(w,T) 

= 2/$| 5sc (q)| 2 coth(^)/ l(g )^, (74) 



Re^,T)= 2j-fJ§ M q) 



2w„ 



ft 2 



(75) 
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It is useful to separate out the zero and finite temperature 
contributions to the mass enhancement by defining 



where 



\ i (T)=\° i +5\ i (T), 



v0 _ dReE?(u,T = 0) 



duj 



UJ = 



SXi(T) 



dReE^(oj,T) 



;=0 



&Re£f(w,T = 0) 



duo 



(76) 
(77) 



(78) 



Here, i £ {1,2} labels the two self-energy contributions 
fEcis l74l75|) . The zero temperature contribution is then 
given by, 



We see that the integrand is negligibly small for phonon 
frequencies il > aT where once again a ~ O(l). In the 
second term, this corresponds to a momentum cutoff of 
q c = a—. Applying these cutoffs we obtain the second 
correction 



5\ 2 (T) = -|Ut 2 , 



where the constant l 2 ~ 0(1) > is given by 



1 

2^2 



dxsech 



(f) 

x I a 2 + x 2 In 



1 

x 2 



(84) 



or 

— ■ 85 

7T Z 



i=l 

2G ri D ^kf 



q dq- 



1 - 



7T fijq) u q 
2 /o(g) r 



(79) 



/o ^<z 

We see that the second term in the round brackets is 
much smaller than 1 for q < Hf- So we neglect this 
term and obtain the zero temperature enhanced mass in 
equilibrium: 

A" « f kF qd q f (q) 



n) Jo 
2G[i 

TT 2 CoV 2 p ' 



(m » r ) 



(80) 



At finite temperature, the effective mass correction gains 
a temperature correction. The first correction is given 
by, 

q 2 dq 



2tt 



C oth ^ - 1 



\2T J 



h(<l)- (81) 



We see that the integrand is negligibly small for q > a— 
where a ~ O(l). Then, employing this momentum cutoff 
we get, 



G (T 

f2 



G 

r 3 r 2 



2tt 

if 

7T Jo 



-KI)- 1 



coth 



(82) 



where ^3 > is a constant of 0(1). The second contribu- 
tion is, 



*«T) = -2 J IfaWlVofa) 



(83) 



In conclusion, the electron mass enhancement at thermal 
equilibrium to lowest order in temperature is 



A(T) 



2G/x 

TT 2 C V 2 F 



GT (T 
-5-/2 I — 



(86) 



2. Out of Equilibrium T = 0, V # 

We begin the noncquilibrium calculation with the sec- 
ond contribution Sf (w, V), 



2 



(87) 



Then the finite voltage correction is given by 



6X 2 (V) = AG I ^fo(q)P 



2tt 2tt' 



n 2 - w 2 



G 



1 p f w 4 ^ 

7T Jo % - 1 



GT 



Here, k 2 ~ 0(1) is a positive constant. For the calcu- 
lation of SXi(V) one needs to go beyond the static ap- 
proximation for the dielectric function as was done for 
the scattering rate calculation. Thus, the correction to 
the effective mass parameter from the first contribution 
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is given by 
6\i(V) = 



Do 2T L T R V 

7r 2 r r 3 

qdqdU,fi{q)^- 



V 



<4 k% {l 



n 2 
W 1 



Do 2T L T R (V 

' 7r 2 r r 2 ^ r 
d g 2 t l t r fvY 



V/c r 2 2 

u q 
qdq-2-2 



(89) 



The electron mass enhancement to lowest order in voltage 
is then, 



X(V) 



2Gn 



-K 2 C V 2 p <", 



GT, 

T 
" 



— fc 2 I - 



(90) 



In conclusion, the electron mass enhancement to lowest 
order in S or -jt, can be summarized concisely by 

X =T : T ^ 0, V = 



2G/i 



GT (X 

X (X) = „ 2 „ „. 2 + 



7r 2 co«|. r;; 



where | x = y | T = > ^ , (91) 
where Kx is a real constant of O(l). 

IV. CONCLUSION 

In this paper we studied a simple model of a 
steady-state electrically-driven two-dimensional electron- 
phonon system. The drive was applied to the metallic 
layer by attaching two 3D leads, which acted both as 
a source of particles (current source) and a heat sink. 
The resultant current was perpendicular to the layer so 
that the heating problem could be avoided. The effec- 
tive theory for the metallic layer was developed using 
the Keldysh path integral method. 

We found that various properties of the electron- 
phonon system is modified in the presence of an out-of- 
plane voltage. Voltage dependences in the Thomas-Fermi 
screening length and in the velocity of the longitudinal 
acoustic phonon mode are presented. The results show 
that both of these quantities can be tuned at will using 
the external voltage. In-plane electron-phonon scatter- 
ing rate and electron mass enhancement were also inves- 
tigated. We showed that electron-phonon scattering can 
be enhanced by voltage at zero temperature. The com- 
puted modification to the electron mass enhancement by 
voltage implies the possibility of tuning the effective mass 
of the electron using voltage. 

The in-plane electron-phonon scattering rate can be 
indirectly measured by observing the voltage-dependent 
piece of the in-plane resistivity, p^ x . The resistivity can 
be measured as a linear response to the in-plane current 
drive. Since the electron-phonon scattering rate occurs 



predominantly in the forward scattering channels, p xx 
gains an additional factor of V 2 compared to the scatter- 
ing rate. Therefore, the finite voltage correction to the 
in-plane resistivity should scale as p xx ~ V 5 . 

Renormalizations to the electron effective mass and 
electron distribution due to external voltage may be ob- 
served in a Shubnikov-de-Haas experiment. In the pres- 
ence of a magnetic field normal to the metallic layer 
(B = Bz), the layer electrons undergo an in-plane cy- 
clotron motion. The electron energies are quantized ac- 
cording to an energy spectrum composed of Landau lev- 
els, which are separated by the cyclotron energy. The 
resultant electron density of states show multiple peaks, 
each peak corresponding to a Landau level. Because the 
separation between a pair of Landau levels depends lin- 
early on magnetic field, the longitudinal in-plane resis- 
tivity oscillates as a function of the magnetic field. This 
oscillation is the well-known Shubnikov-de-Haas effect. 
The difference in the extrema of the oscillations as a func- 
tion of magnetic field decays exponentially as the field is 
decrease d 1 ^ 20 . The cyclotron frequency can be obtained 
by observing the characteristic energy scale with which 
this decay occurs. The effective electron mass, in turn, 
can be found by its direct relationship to the cyclotron 
frequency. 

We also predict that the observed Shubnikov-de-Haas 
oscillations should contain features that specifically re- 
sult from the structure of the zero-temperature nonequi- 
librium electron distribution. Out of equilibrium, the 
distribution function assumes the split-step shape^i, in 
which each step is presumably of different height given 
the asymmetry in the lead-layer couplings. As the mag- 
netic field is varied, both of these steps become resonant 
with a Landau level. As a result, each extremum in the 
resistivity would, in principle, split into two asymmetric 
peaks. 
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APPENDIX A: KELDYSH ACTION FOR PHONONS 



The purpose of this appendix is to begin with the real time action for the Hamiltonian given in section III Al and 
show detailed Keldysh calculations that lead to the effective phonon action Eq|22l from which the inverse retarded, 
advanced and Keldysh phonon propagators can be read off directly. 

Our starting real time action, after the lead electrons have been integrated out, is 



where 



iS[d, d, u,u*] = iS e i[d, d] + iS ph [u,u*] + iS e i- P h[d,d,u,u*], 



iS el [d,d]=i f dtdt'J2d a (k,t')G- 1 (k,t' -t)d a (k,t) 



V ee - 

dt^Y, ~Y d ^ k + 1» ( k ' - <*> t ) d "' ( k '> *)^( k > *)> 

q/0 k,k' 



(Al) 



iS p h[u,u* 
iS e i- p h[d,d,u,u* 



J dtu*(k,t')Do l (k, t'-t)u(k, t), 

-i j dt ^2 Sq<k(k + q,t)d CT (k, i)u(q,t). 



(A2) 



k.q.cr 



In order to carry out the time loop contour integration, we split all fields into two components, forward (+) and 
backward (— ) fields, which reside on the forward and backward branches of the time loop contour, respectively. The 
action then can be written in the form, 

iS [d, d, u, u* , 4>, (j)*] — iS [d+, d+, u + , ul, 4>+, <j)* + ] — iS K [d~, d-,U-,u*_,4>-, </)*_]■ (A3) 

Here, <f> is a Hubbard-Stratonovic field used to decouple the quartic term in Eq lA2l In this basis for the Keldysh 
space (i.e. +,— fields) one obtains four Green functions for every field, one of which can be expressed as a linear 
combination of the other three. Therefore, one often performs a linear transformation of the fields so as to work with 
three independent Green functions. This is known as Keldysh rotation^, in which new fermion fields are defined via 



h = 



h + f- 



and new boson fields are defined via 



b++b_ 



h = 



b q = 



U - f- 



h = 



u - f- 



h = 



U + f- 



, * b* + +b* 
b,,, = 



b* = 

Q 



bl - b* 



2 H 2 cl 2 q 2 

Upon carrying out the Keldysh rotation on the time loop contour action of Eq lA3l we obtain, 



(A4) 



(A5) 



Sf = 2i dt (di(k',t) dl(k',t) 



k,k' 



G^(k,t)S klk [G-Jjf(k,t)« k , k 



Gj 1 (k, t)5k'k 



e<j) cl (k' - k, t) - 5k' -ku cl (k' - k, t) e<\fl (k' - k, t) - <? k ' -kU q (k - k, t) 
&tflQJ -k,t)- <7k'-kU 9 (k' - k, t) e<b cl {k> - k, t) - g k ,- k u cl (k> - k, t) 



4(M) 

4(k,t) 



A W^(k,t) 

h o J { 0«(k,t) 

[^(k*)]- 1 
[flf(k,C 



D^(k,t) 



iSf = 2i / dtJ2((t) cl *(k,t) cf>«*(k,t) 

k 

iS^ = 2i f dtJ2( u cl *(k,t) u q *(k,t) 
J k 

Integrating out the layer electrons in iS^, we get 



(A6) 



u cl (k,t) 
u q (k,t) 



<E> ci (k,w) 
<f> q (k,u) 



(A7) 
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When Eq |A7l is expanded, the total action now becomes, 



ii?fV^<n = 2i£(^*(k,u;) u q *(k,u))( ° , n _! f^'^A 



u cl {k,uj) 
M«(k,w) 



+ 2i£(^*(k,w) 0**(k,w)) 



2tt 



2^£(^ ( k, W ) ^(k,.))f ^ 



^ rf (k,«) 

9 (k,w) 

o n A (k,^) 







k,o; 

+ 2ie^ 5k (^(k ) .) ^(k^))( n , ( ° kjW) 



2i ^| 5k | 2 (^*(k, w ) u«*(k, w ))f nfl(k w) JjV"; 



n A (k,w) 



d (k,w) 
0«(k,w) 

u c/ (k,w) 
w 9 (k,w) 



w d (k,w) 
w) IK u q (k,u) 



(A8) 



The $ fields were defined in Eq[17] and the Keldysh polarization diagrams (II) were defined in Eqs l26l27l Expres- 
sions for the retarded and Keldysh polarization diagrams in the long-time, long-wavelength limit were introduced in 
Eqs l29|30l The coefficients appearing in these approximate expressions are given by, 



ci(V) 
c 2 (V) 

cs(V) 

C A 

Now integrating out the (j> field, we obtain 



r 

E 



Dp \ ■* To M 2 

3tt 



El a 
2T 

a 



DqT 1 

7TM 1 + 



-Do r l 

7T 1 + 



r 2 ' 



- 2i g ( „«.(k,„, ) ( [Df (k ° M)r , w^jji 

o n A (k,w) 



2ie 2 ^| 5k | 2 (^(k^) ^(k^))( nfl(k)a;) nV (k ; 



271 



2 2 n«(k, 



£-e 2 I^(k,c)\ V U A (k,cu) 

-e 2 u K {k,uj) J \n R {k,uj) n*(k,w) 



u cl {k,w) 
tt«(k,u;) 



u c/ (k,w) 
w«(k,w) 



k,w 



2i^( " c/ *(k,w) u q *(k,u) ) 

k,cj 



[^(k,^)]- 1 



e 2 |5k| 







[n- 4 ^)] 2 



i-e 2 n A (k,w) 

[n H (k,^)] 2 n K (k,J) ( i i 
^-e 2 n«(k,^) ? V'Rk^jyp 



- I^k| 2 



/ o n A (k,w) 
Vn fl (k, w ) n K (k,^) 



(A9) 
(A10) 

(All) 
(A12) 



u cl (k,to) 
u q (k,iv) 



(A13) 



We see that our effective phonon Keldysh action obeys the causality structure for bosons. The inverse retarded, 
advanced, and Keldysh propagators can now be read off directly from this action. The results were presented in 
Eqs (2H2S1 
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